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INTRODUCTION 

Hyperfine splitting of tfie ground state in liydrogen was measured already thirty years ago 
up to utmost precision of about 6 x 10^^'^ [||])0- Fot many years this measurement remained 
the most accurate experiment in modern physics. Unfortunately, the theoretical situation 
still remains much less satisfactory in this case, due to uncertainties in the proton-size and 
proton-polarizability corrections. At present, accuracy of the theoretical prediction is on 
the level of several ppm, and further progress is going to be rather slow. This fact limits 
perspectives for testing higher-order QED effects in comparison of theory and experiment 
for the ground-state hyperfine splitting. 

It is possible partly to avoid problems connected with the proton structure, if we study 
the specific difference of the Is and 2s hyperfine structure intervals (Az^i and Az/2), -D21 = 
8Az/2 — Aui. The most accurate measurements of the Is and 2s 0] hyperfine structure 
of ^He"*" yield 

D2i(^He+) = 1 189.979(71) kHz. 



This difference is much less sensitive to the proton structure than Aui and Az/2 separately 
and, therefore, provides better possibilities to test higher-order QED effects by comparison 
with the experiment. 

In the present investigation, we perform a nonperturbative (in the parameter Za) cal- 
culation of the one-loop self-energy correction for the specific difference of the Is and 2s 
hyperfine structure intervals D21 for H-like low-Z atoms. We check consistency of our re- 
sults with perturbative calculations and evaluate the unknown higher-order contribution to 
D21 for hydrogen and He"*". 



I. SELF-ENERGY CORRECTION TO THE HYPERFINE SPLITTING 

In the nonrecoil, point-nucleus limit, the hyperfine splitting of the ns state of a H-like 
atom can be conveniently represented as 
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Au^ = ^[B^ + -D^^\Za)+{-) D^\Za) + ... , (1) 



where Ep is the non-relativistic Fermi energy, i?„ is the relativistic factor, and the functions 
D^'^\Za) represent r-loop radiative corrections. In this work we are interested in the one- 
loop self-energy contribution to the hyperfine splitting, D!^^'^^(Za). This correction can be 
written double expansion in Za and In(Za), 

Dl^^^^^iZa) = aoo + {Za)aio + (^a)'|a22 In'(Za) 
+021 ln(ZQ;) + 020 

+ {Za) [asi HZa) + F^^'iZa)] | , (2) 

where the coefficients aij are listed in Table |, and the function F^^{Za) contains a constant 
entering in order (Za)^ plus all higher-order contributions. 

Now consider the function D^2i'^^ = D^^'^^ — which corresponds to the specific 

difference of the Is and 2s hyperfine structure intervals, 8Ai/2 — Aui. The coefficients aoo, 

('2) SE 

fliO; 022, and 031 are state- independent, and, therefore, the function -D21 can be written 

as 



(3) 



The coefficients 021 and 020 for this difference were investigated in Refs. 
re-evaluated in Ref. (see also references therein), 
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a2i(2s) -a2i(ls) = yln2-7 , 
a2o(2s) - a2o(ls) = -5.221233(3) . 



and recently 

(4) 
(5) 



While the traditional approach to the calculation of QED effects in loosely bound systems 
is based on the perturbative expansion in the parameter Za, it turns out that complexity of 
calculations grows very rapidly in higher orders. As a result, the exact numerical treatment 
of radiative corrections to all orders in Za becomes rather competitive. Despite the fact that 
nonperturbative calculations of the self-energy corrections for light atoms are impeded by 
large numerical cancellations, several accurate results complete in Za have been obtained for 



hydrogen during last few years |]15|-[T8|]. Among them we mention the calculation by Blundell, 
Cheng and Sapirstein |]T5| in which the function Ff^{Za) was determined numerically to all 
orders in Za. For hydrogen, their evaluation yields Ff^{la) = —12.0(2.0). 



II. NON-PERTURBATIVE EVALUATION OF THE SELF-ENERGY 

CORRECTION 



The self-energy correction to the ground-state hyperfine splitting in various H-like ions 
was calculated by a number of authors |I9|- |23|JT5| , |2^ . We note that the ffist two evaluations 
turned out to be not completely correct. In our ffist calculation |T^ we employed the partial- 
wave renormalization procedure. As was found later (see, e.g., |2^), this procedure delivers 
a certain spurious term due to the non-covariant nature of the regularization. This spurious 
term yields an additional contribution of about 3% for bismuth [Z = 83). In the second 
evaluation [^, a certain term in the vertex part was omitted, which contributes on the level 



of about 1% for bismuth (see the related discussion in Ref. ||22|). The self-energy correction to 



the 2s hyperfine splitting was calculated in our previous studies [22,26| for several medium- 



and high-Z ions and, for Z = 83, also in Ref. |2^. In the present investigation, we extend our 



calculations to the low-Z region and increase numerical accuracy, in order to get information 
about terms of order a{Za)^ and higher. 

The self-energy correction to the hyperfine splitting can be conveniently represented as 
a perturbation of the first-order self-energy correction by the hyperfine-splitting interaction 

/X \e\ a. ■ lu X r] , , 

where a denotes the Dirac matrices, and is the nuclear magnetic moment. Taken to 
first order in Vhfs, this perturbation gives rise to the modification of the wave functions, the 
electron propagator, and the binding energy. We refer to the corresponding contributions 
as the irreducible part (Ai^ir), the vertex part, and the reducible part, respectively. The 
irreducible contribution can be expressed as a non-diagonal matrix element of the first-order 
self-energy operator. For its evaluation we use the numerical scheme described in detail in 
Ref. [^. We note that an expression for a wave function modified by the magnetic field 
can be derived analytically both in configuration and in momentum space using virial rela- 
tions for the Dirac equation [^. Corresponding formulas for the Is and 2s wave functions 
modified by the magnetic field can be found in Ref. ||30[| . 

The reducible and the vertex part are conveniently evaluated together. Both of them 
are ultraviolet divergent. In order to isolate the divergences in a co variant way, we separate 
terms in which bound electron propagators are replaced with free propagators. The sum 
of these terms for the reducible and the vertex part is denoted by Ai?^°^ It is evaluated 
in momentum space, where the ultraviolet divergences can be treated in a standard way. 
The remaining part of the reducible and the vertex contribution is referred to as Ai?^]:"^\ 
This part does not contain any ultraviolet divergences and can be evaluated in configuration 
space. We note that both the reducible and the vertex part of AE^l'^'^ possess some infrared 
divergences that cancel each other when considered together. 

Generally, our scheme of the numerical evaluation is rather similar to that used in Ref. 
jTEI . Here we address several features of our calculation which distinguish it from previous 
studies. First, our treatment of infrared divergences in Ai?|J+^ is essentially different from 
that in Ref. . These divergences can be easily separated when bound electron propagators 



are written in the spectral representation, 

As can easily be shown, the infrared divergence arises only when the energy e„ in the spectral 
representation of all internal electron propagators coincides with the energy of the initial 
state Ea- In that case, the integration over the energy of the virtual photon can be carried 
out analytically, and a finite result is achieved explicitly in the sum of the reducible and 
the vertex contribution. So, we divide Ai^^J^^ into two parts. The first one is the finite 
expression derived from the infrared- divergent part, and the next is given by the point-by- 
point difference of the original expression and its divergent part, which is a regular function 
of the virtual-photon energy. 

Some additional complications arise for the self-energy correction to the 2s hyperfine 
splitting, as compared to that for the ground state. For the 2s state, the integration over the 
virtual-photon energy involves the integrand with an additional pole in the first quadrant 
of the complex energy plane. This pole arises from the virtual Is state which is more 
deeply bound than the initial state. When the standard Wick rotation is applied, the whole 
contribution is separated into an integral along the imaginary axis and a pole term. This 
separation leads to a significant numerical cancellation in the low-Z region. An additional 
problem in this case is the presence of a pole due to the 2^3/2 state near the integration 
contour. This pole is separated from the contour only by relativistic effects that nearly 
vanish for very low Z. We avoid both these problems by performing the integration along 
the contour C in the complex energy plane, as shown in Fig. |1]. 

The crucial point of the evaluation, which defines numerical accuracy of the result, is the 
infinite sum over the angular-momentum quantum number of intermediate states, To 
reduce the uncertainty due to the termination of the summation, we extrapolate it to infinity. 
The extrapolation was carried out using two different schemes. First, we use least-squares 
fitting of several last expansion terms to polynomial in Next, we approximate last 

terms of the expansion by where the exponent a is fitted. Convergence of different 

fits as a function of the truncation parameter l^maxl was studied, and the most stable fit was 



chosen. 

Naturally, a reasonable extrapolation of the truncated summation is possible only if 
expansion terms exhibit a clear asymptotic behavior. A problem appears for the irreducible 
contribution, in which the expansion terms change their sign at a certain point \k,\ = \k,o\- It 
is clear that when |ko| comes close to the truncation point |/tmax|; no reliable extrapolation 
is possible. While this problem is present for the ground state if Z is very small, it appears 
much more immensely for the 2s state. We found that by a minor modification of the 
renormalization procedure it is possible to eliminate this unfortunate feature of the partial- 
wave expansion of the irreducible contribution. 

Let us examine the cancellation of ultraviolet divergences in the sum of the zero- and 
the one-potential term for the first-order self-energy correction (for the detailed discussion 
of the renormalization procedure we refer to [Q). The free self-energy operator E*^°^(p) can 
be represented as 

SW(p) ~Sm = -^A,{p,r + Sg)(p) , (8) 

where 6m denotes the mass counterterm, is a divergent constant, p is the 4-momentum. In 
the zero-potential term, the time component of p is fixed at the physical energy of the initial 
state, po = Ea- The renormalized free self-energy operator S^''(p) is finite. Analogously, we 
can write the time component of the free vertex operator as 

T%p,p') = ^A,^^ + TUp,p') , (9) 

where r^(p,p') is finite. Assuming that the operators act on the Dirac wave functions, the 
divergences cancel each other due to the Dirac equation 

(j9,7" - m)Mp) = I (^7Vc(p - p')Mp') ■ (10) 

In the one-potential term, the time components of the 4-momenta p and p' in the vertex 
operator are normally fixed at the physical energy, po = p'o ~ ^a- However, as can be seen 
from Eq. (P), the ultraviolet divergent part of the vertex operator does not depend on them. 
Therefore, we can formally regard e = po = p'q as a free parameter. In this case the one- 
potential term loses its transparent physical meaning but still contains the same ultraviolet 



divergent part. Of course, care should be taken calculating the many-potential term in this 
case. In our approach (see Ref. 123]), the many-potential contribution is calculated as the 



point-by-point difference of the unrenormalized expression and the zero- and one-potential 
terms in configuration space. Therefore, it is sufficient simply to shift the energy in the 
one-potential part from its physical value, Sa —>■ s. 

We found that terms of the partial-wave expansion of the irreducible part do not change 
their sign if the parameter e is chosen as e = {Sa + m)/2. In that case, the expansion terms 
are monotonically decreasing and can be reasonably fitted. A similar procedure can be 
applied to the computation of AE^l~^\ In our implementation, it was used as an additional 
cross-check of the numerical evaluation. 

III. NUMERICAL RESULTS AND DISCUSSION 

The numerical calculation was carried out in the Feynman gauge, for the point nuclear 
model. The partial-wave summation was terminated typically at l^maxl = 40-i-50. The tail of 
the series was estimated by fitting as described above. Various contributions to the functions 

are listed in Tables || and |ni| , respectively. The quoted errors 
originate mainly from the partial-wave extrapolation. The estimated accuracy is essentially 
worse for the 2s state than for the Is state due to a larger contribution of high partial waves. 
For the ground state, our results are in reasonable agreement with the previous calculation 



by Blundell and co-workers [|T^]. However, we note a small systematic deviation of our 
results from those of Ref. |]I5[. For the 2s state, the present evaluation agrees well with our 
previous calculation |2^ of this correction in the medium- and high-Z region. For Z = 49 we 
found D2'^'^^{Za) = —3.555(1), as compared to the point-nucleus result of —3.553(18) p2| . 
The present results for the difference D^2i'^^{Za) show to be consistent with perturbative 
calculations. Fitting our data to the form given by Eq. (|^) with undetermined coefficients, 
we found —3.3(6) for 021 and —5.1(1.0) for 020, which analytically are —3.303 and —5.221, 
respectively. 

Our aim is the higher-order contributions F^^{la) and F^^{2a). In order to get them. 



we have to subtract the known terms of the Za expansion from our numerical resuhs. Of 
course, this subtraction leads to a severe loss of accuracy, especially for very low Z . At the 
present level of numerical precision, the direct determination of the higher-order contribution 
for Z = \ and 2 is not possible. Still, we can evaluate this contribution for higher Z (where 
cancellations are smaller) and then extrapolate it to Z = 1 and 2, as it was done in Ref. 
for the Is state. Aiming this, we tabulate the functions Ff^iZa) and F'^iZa) in a wide 
region of Z starting from Z = 4, as presented in Table We note that for fitting purposes 
the difference F'^iZa) is more appropriate than the function F'^iZa). In this difference a 
number of higher-order contributions cancel each other (e.g., the cubed logarithm of order 
a{Za)'^), and some others have smaller coefficients. 

We start our analysis with the function F2^{Za). The corresponding results are plotted 
in Fig. 1^. We note that the function varies smoothly in a wide range of Z, which again 
indicates that our results are consistent with the first two terms of Eq. (^. Surprisingly 
enough, the general behaviour of the higher-order contribution is rather close to linear in Z, 
even for relatively large nuclear charge numbers of about 30. For very low Z, we observe a 
systematic deviation of our results from the linear behavior, which is, however, well within 
the given error bars. We interpret this as a small systematic error in the partial-wave 
extrapolation of D^2i '^^{^'^) divided by the factor {ZaY when converting to F2'^{Zq). A 
least-squares fit to the numerical results yields F2i{la) = 6.5(8) and F2^{2a) = 6.3(6). We 
estimate the uncertainty of the fitting procedure by studying sensitivity of the result to the 
form of the fit, and assuming that the leading logarithm in the order a{Za)'^ enters with the 



coefficient of about unity. For the introduction into the fitting procedure we refer to [31 



where fitting of nonperturbative results for the one-loop self-energy correction to the Lamb 
shift is discussed in detail. 

Next, we examine the function F^^{Za). The obtained results are plotted in Fig. |^, 
together with the data from Ref. A least-squares fit to our results yields Ff^(la) = 
—14.3(1.1) and Ff^{2a) = —14.5(7). The error of the fit was estimated assuming that the 
leading logarithm of order a(Zaf enters with the coefficient of about unity. Our fitting 
result for Ff^{\a) agrees with the previous calculation of Ref. |jl5| which yields —12.0(2.0). 



IV. SUMMARY 



In the present investigation we evaluated the higher-order contribution to the one-loop 
self-energy correction to the Is and 2s hyperfine splitting of H-like ions. The calculation 
was carried out to all orders in Za. The higher-order correction to the specific difference 
of the hyperfine structure intervals D^x = 8Az/2 — Aui is found to be 0.394(38) kHz for 
3He+, 0.0083(10) kHz for H, and 0.0019(2) kHz for D. For helium, our result should be 
compared with the experimental value D21 = 1 189.979(71) kHz P,^. We see that the higher- 
order contribution evaluated in this work is rather significant for comparison of theory and 
experiment. Still, it should be noted that there is a number of other corrections contributing 
at the same level. For the latest analysis of theoretical predictions for the hyperfine structure 
of hydrogen and He"*" atoms we refer to [jl4[ . 
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TABLE I. Coefficients of the expansion of the function Dn'^^\Za), Eq. (|2|). 



aoo{ns) 
aio(ns) 
022 (ns) 
a2i(ls) 
a2o(ls) 



1/2 



(ln2-f)vr 



-8/3 



(fin: 



37 
36 



17.1227(11) 
-51n2 + ^)7r 



ITIJT 



TABLE II. Individual contributions to Df^{Za). 



z 




aM? 




Df^{Za) 


Ref. |T| 


4 


-0.077320 


2.392195 


-2.073870 


0.241005(25) 


0.24103(1) 


5 


-0.104961 


2.298749 


-2.019760 


0.174028(20) 


0.17405(1) 


6 


-0.134332 


2.210244 


-1.969100 


0.106812(20) 


0.10684(1) 


7 


-0.165112 


2.126160 


-1.921577 


0.039471(20) 


0.03950(1) 


8 


-0.197073 


2.046065 


-1.876933 


-0.027941(20) 


-0.02791(1) 


9 


-0.230047 


1.969610 


-1.834936 


-0.095373(20) 


-0.09535(1) 


10 


-0.263906 


1.896443 


-1.795397 


-0.162860(20) 


-0.16283(1) 


12 


-0.333923 


1.759010 


-1.723003 


-0.297916(15) 




14 


-0.406614 


1.631950 


-1.658572 


-0.433236(15) 




16 


-0.481714 


1.513797 


-1.601149 


-0.569066(15) 




18 


-0.559114 


1.403346 


-1.549961 


-0.705729(15) 




20 


-0.638817 


1.299586 


-1.504357 


-0.843588(15) 


-0.84356(1) 


25 


-0.848878 


1.064059 


-1.411445 


-1.196264(15) 


-1.19621(1) 


30 


-1.077047 


0.854302 


-1.343905 


-1.566650(15) 





TABLE III. Individual contributions to D^^{Za). 



z 






A ■^-,('1 + ') 




4 


-0.06817 


3.95047 


-3.63560 


0.24670(12) 


5 


-0.09171 


3.85376 


-3.58009 


0.18196(10) 


6 


-0.11646 


3.76182 


-3.52821 


0.11714(10) 


7 


-0.14219 


3.67417 


-3.47970 


0.05228(10) 


8 


-0.16871 


3.59045 


-3.43435 


-0.01261(10) 


9 


-0.19588 


3.51034 


-3.39197 


-0.07751(10) 


10 


-0.22365 


3.43355 


-3.35241 


-0.14251(10) 


12 


-0.28070 


3.28902 


-3.28121 


-0.27289(9) 


14 


-0.33956 


3.15520 


-3.21970 


-0.40406(8) 


1 a 
lb 


— U.4UUiO 


o.uouio 


— O.iOlUI 


—U.06o4:l(p) 


18 


-0.46255 


2.91457 


-3.12267 


-0.67065(6) 


20 


-0.52689 


2.80570 


-3.08597 


-0.80716(6) 


25 


-0.69789 


2.56037 


-3.02522 


-1.16274(6) 


30 


-0.88767 


2.34549 


-3.00540 


-1.54757(6) 



TABLE IV. Contributions to Fi^{Za) and Ff^{Za). 



z 


Ff^(Za), this work 


Ff^{Za), Ref. |15| 


F|f(Za) 


4 


-14.71(100) 


-13.85(40) 


7.9(4.8) 


5 


-14.92(43) 


-14.46(21) 


6.7(2.1) 


6 


-15.26(25) 


-14.92(12) 


6.3(1.2) 


7 


-15.51(13) 


-15.283(75) 


5.98(75) 


8 


-15.753(75) 


-15.584(50) 


5.75(50) 


9 


-15.921(49) 


-15.856(35) 


5.60(35) 


10 


-16.188(31) 


-16.113(26) 


5.43(26) 


12 


-16.631(16) 




5.03(13) 


14 


-17.057(10) 




4.712(75) 


Id 


17/1 '7P.(1\ 




4.4oo(^0U j 


18 


-17.896(5) 




4.182(26) 


20 


-18.317(4) 


-18.307(3) 


3.935(19) 


25 


-19.389(2) 


-19.380(2) 


3.336(8) 


30 


-20.524(1) 




2.748(5) 




FIG. 1. The contour C of the integration over the energy of the virtual photon. The poles and 
the brunch cuts of the integrand are shown. 
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FIG. 2. ^21^ 9' function of the nuclear charge number (the upper graph). The same plot with 
the linear part of subtracted (the lower graph). As an illustration to the fitting procedure, 
two examples of fits are shown, the linear and the parabolic one. 




FIG. 3. The function Fi{Z) in different evaluations (the upper graph). The same plot with 
the linear part of Ff^(Z) subtracted (the lower graph). Two examples of different fits are given, 
the linear and the parabolic one. 



